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Abstract. 

Aims. We seek to probe the Galactic bulge IMF starting from microlensing observations. 

Methods. We analyse the recent results of the microlensing campaigns carried out towards the Galactic bulge 
presented by the EROS, MACHO and OGLE collaborations. In particular, we study the duration distribution of 
the events. We assume a power law initial mass function, ^(/i) oc fi~°', and we study the slope a both in the brown 
dwarf and in the main sequence ranges. Moreover, we compare the observed and expected optical depth profiles. 
Results. The values of the mass function slopes are strongly driven by the observed timescales of the microlensing 
events. The analysis of the MACHO data set gives, for the main sequence stars, a — 1.7 ± 0.5, compatible with 
the result we obtain with the EROS and OGLE data sets, and a similar, though less constrained slope for brown 
dwarfs. The lack of short duration events in both EROS and OGLE data sets, on the other hand, only allows 
the determination of an upper limit in this range of masses, making the overall result less robust. The optical 
depth analysis gives a very good agreement between the observed and the expected values, and we show that the 
available data do not allow one to discriminate between different bulge models. 
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1. Introduction 

Gravitational microlensing is an established tool for the 
study and the characterisation of faint compact ob- 
jects located between the observer and the source stars. 
It was originally proposed as a tool for the detection 
of da rk matter in the form of MACHOs (|Paczvnski 
Sea rches towards the Magellanic Clouds by the 
MACHO (lAlcock et all l2000l) and the EROS groups 
( Tisserand et al.l 200?! ) have placed strong constraints on 
the possible contribution of a MACHO po pulation to the 



that of either bulge or disc stars (jGriest et al.l 119911 ). 
so that these studies in principle allow us to constrain 
the inner Galactic structure. In particular, microlens- 
ing observations in this direction have been very impor- 
tant for t he assessment of the Galactic triaxial, bar-like , 
structure (IPaczvnski et"all Il994 IZhao et al.l [l99l Il99d 



Zhao fc Maolll996tlBissantz et al.lll997t ICvudliggg) . 
Recently , the M ACHO (jPopowski et al."2005l). OGLE 



(|Sumi et al I l2006f) and EROS (Hamadachc et all 



dark matter halo (for a discussi on see, e.g.. lMancini et al 



2004 ICalchi Novati et al.ll2006l ). A few results have also 



been obtained with observational campaigns towards M31 
by the POI NT-AGAPE (Calchi Novati et al., ,200a) and 



the MEGA (|de Jong et al.l l2006f) collaborations. On the 



other hand, the Galactic bulge soon pro ved to be an al 



most as interesting target , if not more (|Paczvnskilll99I 



Kiraga &: Paczvnskilll994l ) and indeed, by now, the num- 



ber of observed microlensing events along this line of 
sight is by two orders of magnitude larger than those ob- 
served towards the Magellanic Clouds and M31. In this 
case, any contribution from a dark matter MACHO pop- 
ulation is expected to be extremely small compared to 



collaborations presented the results out of their several- 
year campaigns towards the Galactic bulge. A remark- 
able result is the agreement among the different collabo- 
rations for the optical depth, in accord with theoretical ex- 
pectations tevans fc B elokurovl2002t Bissantz fc Gerhard 
,2002; Han & Gould 200^^ 

While the determination of the optical depth allows the 
study of the density distribution of the lenses, a more de- 
tailed analysis of the shape of the microlensing lightcurves 
carries much information on the parameters of the lenses. 
Of particular relevance is the possibility of studying the 
mass spectrum of the lenses. This approach is based on 
the relationship between the observed event duration, the 
Einstein time, and the mass of the lens oc ^/Jm. Even if 
the exact analytical formula shows a dependence also on 
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other unknown physical parameters such as the distances 
of both lens and source, the relative velocity between them 
as well as the configuration of the particular lens event, 
a few conclusions are made possible by the rather large 
set of observed events at our disposal together with a few 
reasonable assumptions on the space and velocity distribu- 
tions of both lenses and sources. A key step to reach the 
aforementioned agreement of the optical depth between 
theory and experiments has been the acknowledgment of 
the severe blending problem resulting in the choice of re- 
strictin g the sample of s ource stars to the red clump giant 
subset (|PoDOwski|[200lh . In turn, this is essential in the 
framework of a mass spectrum analysis because of the 
bias introduced in the evaluation of the Einstein time for 
blended events. 

The determination of the mass function using the re- 
sults of micro lensing searches has been addressed by sev- 
eral authors. Han fc Gould ( 19961 ) consider a sample of 
MACHO and OGLE events. Through a likelihood anal- 
ysis they determine the slope of a power law mass func 



Jetzer 



tion t o be 2.1 in the mass r ange (0.04 — 10) M0 

(|1994 and lGrenacher et al.l (|l999l ) use the mass moments 
method to place constraints on the lens masses. In par- 
ticular, starting from a sample of 41 MACHO events and 
assuming a Salpeter pr ofile in the mass range (1 — 10) Mq , 
Grenacher et al. ( 1999h constrain the mass function min- 
imum mass and slope below 1 M©, finding 0.012 Mq and 
2.0 respectively. Overall, therefore, there is an agreement 
in attributing a rather large fraction of events to the brown 
dwarf lens population. On the other hand, IPealfJ (|l998h 
finds no compelling evidence for such a contribution, and 
evaluates the slope for a power law mass function in the 
mass range (0.08—2) Mq to be in the range 2.2 — 2.5. All of 
these analyses, we recall, used the complete sample of de- 
tected e vents, not restr i cted t o those with red clump giant 
sources. bissantz et al. ( 20041 ). considering only red clump 
giant sources, find a good agreement with the MACHO 
observed timescale using a mass function with a large 
contribution from the brown dwarf population (with a 
power law slope 2.35 in the mas s range (0.04 — 0.35) Mq) 



lopi 

Wood fc Maol(|2005f ) extend the lZoccah et all (|2000f l slope 
1.3 down to a minimum mass of 0.03 Mq compared with 
the OGLE observed timescale. 

In the present paper our aim is to make use of the most 
recent observational results towards the Galactic bulge to 
study the mass spectrum of the bulge lens population. The 
structure of the paper is as follows. Sect. [2] is devoted to 
the description of the models we use. In Sect. [3] we point 
out a few particular features of the usual microlensing 
quantities upon which we base our method of analysis. In 
Sect. [4] we present and discuss our main results. In Sect. [5] 
we conclude. 

2. Models 

In this section we introduce, describe the features, and fix 
the parameters of our "fiducial" model for the bulge and 
the disc needed to evaluate the microlensing quantities 



that we use in the analysis. Furthermore, we discuss a 
series of changes in the more critical parameters that we 
use to test the robustness of our results. 



2.1. Density distributions 
2.1.1. The bulge 

It is now acknowledged that the Galactic bulge has a 
box-like (tri-axial) structure. In an analysis of clump gi- 
ant stars, IStanek et al.l (|l997h explored several analytical 
distributions to describe the bulge. As a fiducial model 
we use their model E2, which gives the best agreement 
with the observational data, where p(r) = pQexp{—r), 
witlQ r = ^/ix/xoy + [y/yoY + {z/zof, xq = 890 pc 
and axis ratio values a;o : yo : = 10 : 4.3 : 2.8, an 
inclination angle of the bulge major axis with respect 
to the line of sight of a = 23.8° (the bulge is oriented 
with its longer axis pointing t owards us for posi t ive lon - 
gitude values). More recently iRattenburv et al.l ( 2007bl ) 
carried out a similar analysis with a much larger sam- 
ple of stars. As a result, the model E2 is again favoured, 
with axis ratio values suggesting a more prolate structure, 
Xq : yo : zq = 10 : 3.5 : 2.6, and a more restricted range 
of bulge inclination values is given, a ^ (2 4° — 27°). In 



the analysis we keep using the lStanek et al.l (119971) values. 



and we test our results against those of IRattenburv et al 



(|2007bl ). We truncate the bulge at a c orotation radius 
i?c = 3.5 kpc (|Bissantz fc Gerhardll2002l ). The bulge incli- 
nation angle value is still the subject of a somewhat lively 
debate. Values in the range a ~ (10° — 30°) have been 
given by several a uthors, together with different values for 
the axis ratio (e .g. iDw ek et al]ll995l : ISevenster et al.llT999l : 
Picaud fc Robi n 2004), but recently also much larger val 
ues h ave been suggested. Indeed, ICabrera-Lavers et al 



(! 2007t ) (and reference therein) discuss a more complicated 
inner Galactic structure with the co-existence of a dou- 
ble structure, composed of a long 4 kpc) thin and 
lighter bar located at low Galactic latitudes, |&| < 2°, and 
out to high Galactic longitude, with an extremely large 
value for the inclination angle ~ 43°, and a distinct triax- 
ial bulge with smaller incli nation angle, ~ 13°. T he star 
count results of GLIMPSE (jBeniamin et al.ll2005[ ). in the 
I = 10° — 30° range, seem to confirm this result. Such a 
structure may of course give rise to interesting microlens- 
ing signatures, however the currently available data are 
not suitable for its study, as they exclude the Galactic 
plane region and are mostly restricted to events observed 
at small Galactic longitudes. (The EROS collaboration 



()Tisserand et al.l 120071 ) evaluates a bulge orientation an- 



gle of 49° ± 8° even if they observe only the region out to 
|Z| ^ 10° and do not observe any field for \b\ < 1° and only 
a very few in the band out to |&| ~ 2°.) T he issue of the 
bulge inclination has also been discussed in IWoodI (|2007l ) 



^ Here, as in Sect. l2.2l where we discuss the velocity distribu- 
tion, the coordinates x, y, z indicate the principal axes of the 
component considered, namely, either of the bulge or of the 
disc. 
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in the framework of an analysis of the micro lensing optical 
depth As a n alternative bulge distribution. iHan fc Gould 
( 19951 2003f) use the model G2, fa voured by an anal ysis of 
the COBE DIRBE observations (|Dwek et alJIlQQSl ). that 
they cor rect fo r smal l (r < 700 pc) galactocentric distance 
with the KentI ( 19921 ) model. We compare these two mod- 
els to the observed optical depth. The total bulge mass is 
usually evaluated in the range MhnVs ^ (1 — 2) 10-^°Mf ;^ 
( Blumlll995l : IZhao et all Il996l : IPehnen fc BinnevI [19981 ). 
Lacking any compelling constraints we choose to normalise 
the bulge distribution to the observed value of the mi- 
crolensing optical depth (Sect. 14. ip . Throughout the pa- 
per we use Rq = 8 kpc as the value for the distance to 
the Galactic centre. 



2.1.2. The disc 

The profile of the disc distribution is better constrained 
than that of the bulge, although the value of the param- 
eters that characterise it is subject to de bate. In order 
to par ametrise the model we follow closely iHan fc Gould 
( 20031 ) who use a sech^ (exponential) profile for the thin 
(thick) components and normalise the distribution so as 
to obtain a local s tellar density of Sn = 36 Mq pc~^. We 
note however that Han fc Gouldl (2003) attribute a rather 
large density fraction to the "thick" disc, whereas this 
component is usually reported to contr ibute only to a mi- 
nor f r action of the overall density (e.g. iDehnen fc Binnev 
19981 Vallenari et al.llioO^ . Fo r our fiducial model w e as- 
sume, as compared to that of Han fc Gouldl (2003), the 
extreme case where we set to zero the thick disc con- 
tribution (in their notation, we use P — instead of 
/? = 0.565); moreover, we fix the value of the local disc 
density in agreement with the normalisation of the disc 
mass function (Sect. 12. 3[). We have then t ested our 
suits using th e value s of iian fc Gouldl (|2003h and also the 
iFreudenreichI ( 19981 ) profile characterised by a decreasing 
density towards the Galact ic centre. As already pointed 
out by iHan fc Gouldl HoO^), we find that our results do 
not depend significantly upon the disc model as the bulge 
component gives by far the dominant contribution to the 
observed events. 



2.2. Kinematic models 



To evaluate the microlensing rate we have to specify the 
velocities of the components involved. For bulge and disc 
stars we take into account both the bulk and random 
components of motion. For the former, for both disc and 
bulge we assume a solid body rotation out to i?cut and 
at outer radii a flat rotation with i?cut ~ 2 (1) kpc 
and Knax = 220 (50) km/s for the disc (bulge) com- 
ponent r espectively. As the bulge value is less well con- 
strained (lBlumlll99£tlDehnenll2000l:|Bissantz et al.l l2003l; 



Minchev et al.ir2007tfRich' et al.ll2007n . we have tested our 
results varying the bulge component by 30% to larger and 



smaller values. Furthermore, we take into account the so- 
lar motion. 

For the random component, we assume the velocity 
distributions to follow an anisotropic Gaussian profile. For 
the disc dispersion we use ax = 20 km/s and we consider 
a linear increase towards the Galactic centre for the re- 
maining components with {ay,az) = (30,20) km/s and 
{cryiC^z) = (75,50) km/s in the loca l neighbourhood an d 
at the Galactic centre respectively ( Han fc Gouldl [l995 ) . 
For the bulge, whose velocity dispersions are not as well 
constrained, we consider two somewhat opposite c ases. 
As a first approach, we follow iHan fc Gouldl (|l99,5[ ) and 
fix the dispersion values using the virial t heorem as ap- 
plied to the bulge distribution (|Bluml Il995h . For our fidu- 



cial model we obtain cr^.y^z — (112.5,86.1,72.1) km/s. In 



Sect. 14.2) we investigate the effects of changes with respect 
to our fiducial model. Whenever we modify either the cen- 
tral density or the pattern speed of the bulge component, 
we adjust the dispersion values of the bulge according to 
the prescription of the virial theorem. In only one case 
we arrive at rather significant differences (beyond a few 
pe rcent), namely , when we consider as a disc model that 
of iHan fc Gouid (|2003h . Indeed, in that case, our evalua- 



tion of the bulge central density decreases by about 25% 
(Sect, mi]), implying (7x,y,z = (96.3,74.3,62.8) km/s. 

As a second estimate, we make use of recent obser- 
vation al results ( Kozlowski et al. 1 120061: iRattenbury et al 
2007al ) an d use (ai,ah) = 3.0, 2.5 ma s yr~^, with ctj, 



110 km/s (jBinnev fc Merrifieldl IT998h . For a bulge incli- 
nation of a = 23.8° and Rq = 8 kpc we then evalu- 
ate the dispersion along the bulge principal axes to be 
(109.4,114.8,94.8) km/s. 



re- 2.3. IVIass functions 



The main aim of the present work is to analyse the mi- 
crolensing events to place constraints on the mass function 
of the bulge stars. Zoccali et al.l ( 2000l ) study the bulge 
mass function in the range (0.15 — 1) Mq, finding a good 
fit to the data with a IMF power law, ^(/i) cx with 
a = 1.3±0.1. They also propose a power law with a change 
of slope at 0.5 Mq and a ~ 1.4, 2.0 respectively below 
and above this threshold. Overall, thi s result is comp atible 
with the previous analysis of Holtzman et al.l ( 19981 ). The 
more difficult part of the mass spectrum to be explored is 
the low mass tail, including very low mass main sequence 
stars and the brown dwarf range. In their analysis of mi - 
crolensing e vents towards the bu lge, iHan fc Gouldl (|2003f ) 



extend the Zoccali et al. ( 2000l ) mass function down to 
well below the hydrogen mass burning limi t, at 0.03 Mp). 
and t h e same was d one more recently by IWood fc Mao 
(l2005l ). lGouldl ('2000') describes how to treat remnants, as- 
suming that all of the stars with mass above 1 Mq have 
by now entered the rem nant phase. G iven a slope of the 
IMF in this mass range, iGouldl ^Qod) proposes a = 2, it 



is then possible to evaluate the number and mass fractions 
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due to each of these components (white dwarfs, neutron 
stars and black holes). 

Following the previous results, we assume a power 
law mass function for both the brown dwarf and the 
main sequence ranges. We introduce two parameters, the 
slopes aBD, ckMS in the mass ranges (0.01 — 0.08) M0, 
(0.08 — 1.0) Mq respectively, that we want to constrain. 
According to analyses carried out for the disc, the slope 
should chang e below the hydrogen burning limit (e.g. 
Kroupal2007 ). We also test the effects on our results of two 
changes on the bulge mass function, namely we introduce 
a slope change at .5 Mq, using a — 2 above this limit 
( Zoccali et al. I l2000l) . and we move the lower brown dwarf 
limit from 0.01 to 0.04 Mq. We follow |Gould| ^qM) to 
deal with the remnants, with (0.6,1.35,5.0) Mq taken as 
the mass values for white dwarfs, neutron stars and black 
holes respectively. Besides the value aicm = 2, in this mass 
range we will further test our result with the higher value 
arem = 2.7, as suggested by disc results. Note that for ev- 
ery pair of values {ubd, ckMs)i the number and the mass 
fractions of the various lens components change accord- 
ingly. 



Fo r the disc mass function we closely follow iKroupa 
(|2002l l2007l ). with a power law with slopes 0.3, 1.3, 2.3 
in the mass ranges ((0.01 - 0.08), (0.08 - 0.5), (0.5 - 
1.0)) Mq (the low value in the br own dwarf region is 
in agreement with I Allen et al.ll2005l) . and normalisation 
Io 687 '^l^ ~ stars pc~^. To account for the 

re mnant contribut ions we use the density values reported 
in IChabrieil (|2003f) . to obtain E''™^ = 3 Mq pc'^. This 
fixes the overall local density for our fiducial disc model 
to 4.4 10^ Mq kpc~^ 

3. Analysis: the microlensing quantities 

Our main tool of investigation is the rate of microlens- 
ing events F, that carries the information of the num- 
ber of events per time interval, whereas the microlens- 
ing optical depth, r, is the instantaneous probability 
of a star being magnified above a given threshold (e.g. 
Roulet fc Mollerach 19971 ). Through the analysis of the 
differential rate, given the efhcicncy of the experiment, 
one can analyse the distribution of the relevant microlens- 
ing parameters as well as evaluate the number of expected 
microlensing events. 

The microlensing rate (|De Ruiula et al. 1 119911 : iGriestI 



199lh depends, for both sources and lenses, on the density 
and velocity distributions, on the lens mass function and 
on the microlensing configuration. Once the theoretical ex- 
pression for the differential rate is obtained fAppendix E]) . 
to compare with the results of a given experiment, we still 
need to specify the efficiency of the analysis, usually pro- 
vided as a function of the microlensing timescale together 
with the value of the maximum impact parameter allowed. 

Throughout the paper we will only consider the sim- 
pler microlensing event configuration, point-mass lens 
and source with uniform relative motion between lens 



19861 ) . The effects of non-standard configuration events for 
the evaluation of the microlensi ng quantitie s have been 
the object of a detailed study bv lGlicensteinI (j2003l ). The 
largest changes are to be expected for binary caustic cross- 
ing events, but these represent only a very small fraction 
of the overall set so that the modifications in the evaluated 
quantities should not exceed 10 %. 

4. Results 

4.1. The optical depth profiles 

The agreement among the different collaborations 
(MACHO, EROS and OGLE) on the value of the ob- 
served optical depth, and its agreement with theoretical 
models is, as already noted, a significant result of the mi- 
crolensing searches towards the Galactic bulge. We take 
advantage of this result by making the choice to normalise 
the bulge central density to the observed value of the 
optical depth. As a fiducial value we take the result re- 
ported by the MACHO collaboration tow ards the "CGR" 
( "Cen tral Galactic Region" , defined in iPopowski et al. 



(j2005[ ) as 9 out of the 94 observed fields nearest to the 



Galactic centre), namely r = 2.17l[;;3^ 10 "^ for (lb) 



1°. 50, —2°. 68. For our fiducial model, the Stanek et al 



(|l997l ) model E2, this gives us a central bulge density of 
po = 9.6 10^ Mq kpc"'^, corresponding to a bulge mass 
out to 2.5 kpc of 1.5 10^° Mq (this is strictly the mass due 
to possible lenses). For the model G2 we obtain instead 
Po = 2.4 10^ Mq kpc"^ and a mass of 1.4 10^° Mq. 

Having normalised our model to the optical depth ob- 
served along a given line of sight, next we have to test the 
optical depth profile against the observed one, given that 
the observed events are spread over ^ 4° in Galactic lati- 
tude and ~ 10° in Galactic longitude. Indeed, the optical 
depth profile depends strongly on the line of sight, in par- 
ticula r on the Galactic latitude (e.g. lEvans fc Belokurov 



2002t ). In Fig. [T] we show the optical depth profile for the 



models E2 and G2. We note the larger gradient along the 
Galactic latitude for the first model. 

To gain fu rther insight on their results, the EROS 



20061) 



studied the rela- 
giving the empirical expression r — 



and source, the so-called Paczynski lightcurve (Paczvnski 



collaboration (jHamadache et al. 
tion T — rib) 
N exp[— a(|6| — 3°)], where a, N are to be determined 
from the observa t ional data. For the EROS data set, 
iHamadache et all (|2006l ) find N = 1.62 ± 0.23, a = 
0.43 ± 0.16. As a theoretical prediction, given the 
EROS observational setup, for the E2 (G2) models we 
find N,a = 1.77,0.52 (1.87,0.37) respectively. If we carry 
out the same exercise considering either the MACHO or 
the OGLE observational setup we find the values N,a = 
1.60, 0.56 (1.72, 0.39) and N,a = 1.81, 0.51 (1.94, 0.34) re- 
spectively. Overall, the E2 and the G2 model predictions 
are both consistent with the observed values. 

The previous analysis has to be carried out taking 
bins in the Galactic latitude, averaging over the Galactic 
longitude for the observed fields. This way, however, one 
misses the information of the (albeit smoother) variation 
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Fig. 1. Optical depth profiles as a function of the Galactic 
coordinates l,b for the bulge models E2 (top) and G2. 
The contours of the 9 MACHO CGR fields are shown. 
The optical depth is normalised to the value of the 
observed r at position l,b = 1°. 50,— 2°. 68 (see text 
for details). The profiles are drawn at values of t = 
(0.3, 1. 0, 1.5,2.17,3.0) IQ-^ Overall, the 94 MACHO 
fields (jPopowski et al.l 120051 ) extend in the ranges 0° , 8° 
and —2°, —10° in Galactic longitude and latitude r espec- 
tively. The 66 EROS fields (|Hamadache et al.ll2006[ ) cover 
two regions at both positive and negative Galactic lati- 
tude, / ~ (8°, -6°), b ~ (-2°, -6°) and I ^ (6°, -4°) , 
b - (2°, 6°). The 30 OGLE fields analysed in ISumi et al 



(200a) cover a smaller region near the Galactic centre 
spreading only slightly beyond the MACHO CGR fields. 



of the optical depth profile along the Galactic longitude. 
Moreover, of course, we are comparing the expected op- 
tical depth to the EROS observed values only. As a dif- 
ferent approach, we propose to take bins, instead, in the 
expected optical depth, to be compared with the observed 
one as evaluated for each observational campaign. To per- 
form this analysis we first need to evaluate the observed 
value of the optical depth, and therefore the number of 
sources stars, in each chosen bin, whereas this number is 
known per field. As a first order approximation we con- 
sider the number of source stars in a given fraction of a 
field to be proportional to its area. A bin in the theoreti- 
cal optical depth delimits a region in the Galactic plane. 
We choose the bin sizes so to get a (roughly) equal num- 
ber of observed events in each bin. For EROS (MACHO) 
we tried with both 5 and 10 (3 and 5) bins, resulting in 
very similar results; for OGLE we use 3 bins. In Fig. [2] we 
show the observed optical depth as a function of the ex- 
pected one, for the EROS and MACHO and OGLE data 
sets (with 5, 3 and 3 bins respectively) and both models 
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Fig. 2. The observed and the expected optical depth for 
the bulge models E2 (top) and G2 (see text for details). 
EROS, MACHO and OGLE data are the empty, filled 
circles and stars respectively. The solid line is the y = x 
line. 



E2 and G2. We find very good agreement for both mod- 
els with the three data sets. Indeed, if we fit the relation 
Tabs = o, ■ Tth, considering the 11 points of the three data 
sets together, we get a — 0.9±0.1, for a reduced = 1 
for both models. 

A possible way to disentangle the different models 
would come from an independent normalisation of the 
bulge mass. Indeed, according to our choice, the expected 
optical depth is made to coincide at the CGR MACHO lo- 
cation, therefore, even if different, the two profiles remain 
rather near each other along the observed fields. The other 
option would be to observe events over a larger area of the 
sky. The more interesting region to be explored being that 
closer to the Galactic plane. 

4.2. The Galactic Bulge IMF 

The microlensing rate, as discussed in Sect. [3l is an ef- 
ficient tool for the analysis of the characteristics of the 
microlensing events. Here we focus on an analysis of the 
timescales provided by the current observations. Indeed, 
as outlined in the Introduction, though degenerate with 
other unobservable quantities (distances and relative ve- 
locity between sources and lenses) the dependence on the 
lens mass makes the timescale a valuable source of in- 
formation on the mass function of the lens population. 
Due the above-mentioned degeneracy, one needs a rather 
large number of observed events to deal with them statis- 
tically. The current observational results begin to provide 
such data set, 62 events from the MACHO collaboration 
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(IPopowski et al.ll2005l). 120 events from the EROS collab- 
oration ( Hamadache et aP 2006) and the 32 events from 



the OGLE collaboration (jSumi et al.ll2006l ). Note that we 



choose not to consider the different data sets together, 
rather, we carry out independent analyses and then com- 
pare the results. 

The model, as described in the previous section, to- 
gether with the microlensing event geometry and the ex- 
perimental apparatus, summarised in the reported detec- 
tion efhciency usually given as a function of the duration, 
£ = f (^e), provides us with the expected number den- 
sity of the microlensing events. Allowing for the Poisson 
nature of th e process we can write down the likelihood 
(|rxouldll2003h . as a function of the free parameters of our 
model, 



L(aBD, aMs) = exp(-iVoxp) 



dr 



i.£ 



(1) 



Here A'cxp is the overall expected number of events, to be 
evaluated by integrating out the differential rate taking 
into account, besides the detection efficiency, the number 
of sources and the overall duration of the experiment. In 
particular it results A'oxp = -^exp(aBD, aivis)- 

As outlined in Sect. 12.31 we take as free parameters the 
slopes of the IMF in the brown dwarfs and main sequence 
ranges, aBD, cimSi that we want to estimate. To evaluate 
the likelihood, we sum the disc and the bulge contribu- 
tions, and for each the contribution of the brown dwarfs, 
main sequence and remnants lens populations. 

Finally, to estimate the confidence levels, we evaluate 
the probability distribution P (aeD, c^ms) by Bayesian in- 
version using a flat prior on both the parameters. 

It is useful, for our purposes, to take the sample of 
MACHO CGR events as a "fiducial" sample. This pro- 
vides us with a more homogenous, but still quite large, 
set of events all located in a region small enough to 
make any possible spatial dependence, that we may not 
have correctly reproduced within our model, almost irrel- 
evant. Furthermore, the CGR allows a more straightfor- 
ward comparison among t he different data sets. 

The 66 EROS fields (jHamadache et all l2006l) cover 
a rather larger region in the plane of the sky than the 
MACHO fields, both at positive and negative Galactic 
latitude. For comparison with the MACHO CGR sample, 
we select the 5 fields (5,8,607,610,611) whose location is 
roughly coincident with that of the MACHO CGR fields, 
and where 18 of the 120 events reported by the EROS 
collaboration are located. 

Finally, we observe that the l ocation of the 20 fields 
used by OGLE in their analysis ( Sumi et al. 20061 ) only 
slighly exceeds the CGR. 

For the observed distributions for both the MACHO 
and the EROS data sets, there is an increase in dura- 



^ This is the so-called "extended maxim um likelihood" , first 
proposed by Fermi (for a discussion see e.g. iBarlowll 19891 ). that 
is appropriate in experiments where the number of events is 
itself a random variable. 



Table 1. The average observed duration, (ts) (days), 
of the microlens i ng ca ndidates r eported by the MACHO 
(|Popowski et al.l l2005f). ERO S (jHamadache et all l2006t) 
and OGLE ( Sumi et al. 20061) collaborations. In the first 
row, for MACHO and EROS we report the result within 
the CGR (see text for details). For each data set, in the 
right column we report the average weighted by the in- 
verse efficiency. 





MACHO 


EROS 


OGLE 


CGR 


19.5 15.0 


25.9 22.4 




all set 


28.0 20.0 


32.9 28.3 


32.8 28.1 



tion moving from the smaller sample in the CGR to the 
complete data set. In Table [T] we report the average ob- 
served durations, both uncorrected for the efficiency and 
weighted by the inverse efficiency, the latter quantity al- 
lowing a more straightforward comparison between the 
different data sets. 

4.2.1. The analysis within the CGR 

The main result of the present paper is shown in Fig. [3l 
From the maximum likelihood analysis we show the con- 
tours of equal probability in the q;bd, Q^ms parameter 
space. Here we consider the sample of the 42 MACHO 
events observed within the CGR. 

The data better constrain the IMF slope in the main 
sequence range than in the brown dwarf range. As for 
the IMF parameters, at maximum probability we get the 
values q:bd = 1-6, aivis = 1-7. The corresponding bulge 
mass fractions are ~ (21%, 56%, 17%, 4%, 3%) for brown 
dwarfs, main sequence, white dwarfs, neutron stars and 
black holes, respectively, for an average mass of 0.1 Mq. 
Note the rather high brown dwarf fraction, indeed within 
the 34% level it does not decrease below ~ 20%. Overall 
the bulge contributes about 80% of the events (this result 
confirming the statement made about the only relative 
importance of the disc contribution, Sect. I2.1.2P and the 
event fractions due to the different lens populations are 
- (29%, 57%, 11%, 2%, 1%). In Fig. H we show the one 
dimensional probability proffie P {aBu) and P(aMs)- As 
already mentioned, the Q!ms distribution turns out to be 
more peaked, with q;ms = 1-7 ± 0.5 and q;bd = 1-6 ± 1.0. 

In Fig. [3] we show the lines of equal value of the ex- 
pected duration t^ superimposed on the likelihood prob- 
ability contours. As it is apparent from the plot, the lines 
of degeneracy in the parameter space ajvis ~ Q^bd that are 
found in the probability contours are driven by the du- 
ration (within the innermost 34% probability contour the 
dispersion of the expected duration is only about 5%). In 
particular, we observe that expected shorter durations are 
associated with steeper mass function. This is expected, 
of course, because of the relationship between the dura- 
tion and the mass of the lens. This correlation is relevant 
in order to properly understand the variations we find in 
the evaluated slopes of the mass function for either sets of 
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Fig. 3. Probability isocontours with 34%, 68% and 90% 
regions in the aBD, aivis plane. ofBD, aus are the slopes 
of the power law IMF of the Galactic bulge lenses, in the 
brown dwarf and main sequence range, respectively. The 
dashed lines are the lines of equal average expected event 
durations, for the values of 10, 15, 20, 25 and 30 days. 
Larger values of the duration are found for smaller values 
of the IMF slopes. Here the set of 42 events reported by 
the MACHO collaboration in the CGR is considered. 

data with different duration distributions or for different 
models. 

We now compare the results we obtai n using the sam- 
ple o f MACHO microlensing candidates ('Popowski et al 



20051) with that of the ER OS (.Hamadache et al..,.2006i) and 
of the OGLE (ISumi et al.ll20oir collaborations Following 
the previous discussion, we first consider the samples re- 
stricted to the inner Galactic region, namely we use 18 
out of the 120 EROS microlensing candidates and the full 
set of the 32 OGLE microlensing candidates. 

As shown in Fig. [5] (top panel), the analysis over the 
EROS data set allows us to determine the maximum for 
the IMF slope in the main sequence region, roughly consis- 
tent with that found using the MACHO data set, but does 
not reveal any lower limit in the brown dwarf range. This 
arises because of the different distribution of the observed 
timescale. In particular, the exp lanation may be traced 
back to the lack (already noted in Hamadache et al1l2006f ) 
of very short duration events, say below 5 days, within the 
EROS data set (both in the restricted sample of 18 events 
we consider here and in the full data set). While this differ- 
ence docs not significantly affect the results on the optical 
depth, in the present analysis this turns out to be very 
relevant. The analysis performed on the OGLE data set 
provides a qualitatively similar result. In agreement with 
the previous discussion, we note that the somewhat lower 
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Fig. 4. Probability distribution for the Galactic bulge 
power law IMF slopes in the brown dwarf (top) and main 
sequence mass ranges. Here the set of 42 events reported 
by the MACHO collaboration in the CGR is considered. 
The probabihty fines of 16%, 50% and 84% are indicated. 



value for aMS is a consequence of the higher average ob- 
served timescale. 

The above analysis clearly shows the extent to which 
short duration events are essential to constrain the lower 
tail of the IMF. We further address this question, also 
to compare the results we obtain with the different data 
sets, by mean of the following analysis. We set to zero the 
efficiency below a given threshold, in particular £{t-E < 
5 d) = and at the same time we exclude from the analysis 
those observed events with t^ < 5 d, namely, the 6 events 
from the MACHO sample. 

The likelihood contours we obtain for the MACHO 
data set are shown in the bottom panel of Fig. [5l 
Comparing with Fig. 0] we see that, as for the EROS and 
the OGLE data sets, the brown dwarf slope is no longer 
bounded at its lower end, while the main sequence one 
peaks roughly in the same region. Carrying out this anal- 
ysis for the EROS data set we find an almost identical 
result to that shown in the top panel of Fig. O while for 
OGLE we find a somewhat different behaviour. In that 
case a lower bound for Q!bd appears, at least for the in- 
nermost 34% contour, but at the same time the contours 
become unbounded at the upper end. A similiar behaviour 
is also observed when we move the lower limit on the lens 
mass from 0.01 to 0.04 M0, Sect. I4.2.3i and this can be 
understood, as we are excluding the duration range where 
the microlensing rate of very large brown dwarf slopes 
peak. 
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Fig. 5. Probability isocontours with 34%, 68% and 90% 
regions in the aBD, aivis plane. asD, oiMS are the slopes 
of the power law IMF of the Galactic bulge lenses, in the 
brown dwarf and main sequence ranges, respectively. From 
top to bottom, the results of the analysis for the EROS and 
OGLE data and the results of an analysis of the MACHO 
data where we set the efficiency below is < 5 d to zero (see 
text for details). For the MACHO and EROS data sets we 
restrict the analysis to the subset of events observed in 
the inner Galactic region. 



4.2.2. The analysis of the complete data set 

The analysis of the complete data sets confirm our pre- 
vious conclusions. In Fig. [6] we show the probability con- 
tours for both the full set of events of MACHO and EROS 
(for OGLE, the results obtained with the full data set are 
shown in Fig. [5]). For the MACHO data set we evaluate the 
slope in the main sequence range to be q;ms — 1.6 ±0.4, in 
agreement with the previous result. With respect to Fig. [3] 
and Fig. [5] we observe, however, the maximum likelihood 
contours moving towards somewhat smaller values of the 
IMF slopes. This is of course to be attributed to the in- 
crease in the observed duration (Table [1]). Comparing to 
Fig. [3 for both MACHO and EROS data sets we observe 
a shrinking in the probability contours due to the much 
larger sample of events used in the present analysis. 

4.2.3. The IMF: a test against the fiducial model 

We now focus on the possible systematic effects resulting 
from a change in the characteristics of our fiducial model 
(Sect. [2]). We carry out this analysis using the MACHO 
data set only within the CGR. To give an indication of 
the goodness of the model we use a Kolmogorov-Smirnov 
(KS) test comparing the expected to the observed dura- 




0.5 1 1.5 2 2.5 3 3.5 4 

MACHO data ^ms 




1.5 2 2.5 

EROS data 



a 



MS 



Fig. 6. Probability isocontours with 34%, 68% and 90% 
regions in the q;bd, c^vis plane. Q!bd, Q^ms are the slopes 
of the power law IMF of the Galactic bulge lenses, in the 
brown dwarf and main sequence ranges, respectively. The 
full set of events for MACHO (top) and EROS data sets 
are considered. 



tion distribution (we report its significance level, fcs, as a 
disproof of the null hypothesis that the distribution are 
the same, such that a low value of ks indicates a poor 
agreement between the expected and the observed distri- 
bution). Such an analysis makes sense because the varia- 
tions in the expected timescale do not exceed ~ 5% across 
the CGR. This allows us to carry out the KS test by eval- 
uating an average rate, summing the rate observed to- 
wards the different fields with a weight given by the num- 
ber of source stars that we compare to the observed dura- 
tion distribution. As test models we consider the following 
(Tabled]): model 2 : we change the b ulge dispersion ve- 
locity to the iRattenburv et al T (|2007a[ ) values (Sect.d^]); 
model 3-4 : we change respectively downward and up- 
ward the bulk rotation velocity of the bulge (Sect. 12. 2[) : 
model 5 : we change the parameters of the disc den- 
sity profile according to lHan fc Gouldl (|2003f) (Sect.[lIII2D; 
model 6 : we chang e the disc density profile according to 
iFreudenreich ( 19981) (Sect. \TT7^: model 7 : we change the 



bulge scale lengths according to IRattenburv et al.l (|2007bD 
fSect. [^TTTT]) : model 8 : we change the bulge IMF slope to 
arom — 2.7 (Sect. 1131); model 9 : we change the bulge mass 



function introducing a second slope in the main sequence 
range (Sect. 12. 3^ : 10) we change the mass lower limit in 
the brown dwarf range to 0.04 M© (Sect. [231). In Table [2] 
we report the results: for each model we give the evaluated 
«MS parameter out of the P(aMs) distribution (with the 
16%, 50% and 84% bound) and the KS significance level. 
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Table 2. Results of t he maximum likelihoo d analysis on 
the MACHO data set (|Popowski et al.ll2005l) for our fidu- 
cial model and for the different models discussed (Sect.[2|). 
ajvfs is the slope of the power law IMF in the main se- 
quence range, ks is the Kolmogorov-Smirnov significance 
level for the null hypothesis that the expected and ob- 
served distributions are the same. 



model 




QMS 




ks 


model change 




16% 


50% 


84% 






1 


1.18 


1.74 


2.20 


0.39 


fiducial model 


2 


0.86 


1.40 


1.80 


0.29 


bulge velocity 


3 


1.16 


1.73 


2.19 


0.39 


bulge velocity 


4 


1.17 


1.74 


2.20 


0.39 


bulge velocity 


5 


1.32 


2.10 


2.90 


0.41 


disc model 


6 


1.13 


1.67 


2.10 


0.40 


disc model 


7 


0.97 


1.51 


1.92 


0.33 


bulge model 


8 


0.95 


1.56 


2.06 


0.40 


mass function 


9 


0.90 


1.62 


2.21 


0.40 


mass function 


10 


1.86 


2.19 


2.60 


0.26 


mass function 



For the different models the likelihood maximum 
moves on the aBD-aMS plane so as to always peak around 
the same expected timescale, with the resulting mass func- 
tion slopes changing accordingly. The largest variation 
downward, aMS ~ 1.4, is found for model 2 as an ef- 
fect of the increased bulge velocity dispersions. Note the 
large value we obtain for model 5, we find oms ~ 2.1. 
Here two different effects push in the same direction to- 
wards a steeper mass function, namely a smaller bulge 
contribution and a decrease in the bulge dispersion veloc- 
ity. The qualitative shape of the likelihood contours does 
not change for any of the models except the last. Here, as 
an effect of the increase of the minimum mass value in the 
brown dwarf range, from 0.01 to 0.04 M0, the probability 
distribution for aeo becomes unbounded at its upper end. 
Correspondingly, we also find a steeper mass function and 
lower KS significance level. Overall, the variations we find 
for aivis for the different models we have tested do not 
exceed the statistic uncertainty we have in our fiducial 
configuration. This is in agreement with the KS analysis, 
according to which we obtain acceptable results for all the 
models we consider. 

In Fig. [7] we show the cumulative distribution for the 
sample of the 42 MACHO CGR events together with the 
theoretical cumulative distributions for the fiducial model 
and models 2 and 10, for which we obtain the small- 
est and the largest values for the main sequence slope 
(q^ms = 1-4, 2.2, respectively) and the worst agreement 
according to the KS test. Besides the lack of observed 
events at ^ 20 d, we note in particular the very good 
agreement with short duration events for both the fidu- 
cial model and model 2 and the better agreement with 
long duration events for model 10. 
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Fig. 7. Cumulative duration distributions (P) for the ob- 
served events, together with the the theoretical predic- 
tion for three different models. Solid, dashed and dotted 
lines are for the fiducial model and models 2 and 10, re- 
spectively (see text for details) . Here the set of 42 events 
reported by the MACHO collaboration in the CGR is con- 
sidered. 



4.3. The expected number of microlensing events 

Besides the study of the duration distribution, the anal- 
ysis of the microlensing rate allows one to evaluate the 
number of expected events. Through the analysis we have 
normalised the bulge central density, once having fixed 
that of the disc, by using the observed value of the opti- 
cal depth. Because of the relationship between the optical 
depth and the microlensing rate, through the event dura- 
tion, we may therefore expect to find a good agreement 
between the observed and the expected number of mi- 
crolensing events. Indeed, even if the number of expected 
events varies by almost a factor of 3 across the asD — ctMS 
parameter space we explore, we find a fair agreement. For 
the MACHO data set our prediction is compatible within 
1 (7 to the observed value; we find 38 and 54 events com- 
pared with 42 and 62 events, in the CGR and the complete 
data set, respectively. For the EROS and OGLE data sets 
we arrive at an even better agreement, with an expected 
number of 118 (31) compared with 120 (32). These figures 
do not vary significantly (at most by ~ 2 events) within 
the innermost 34% probability contour. 

4.4. The blending issue 

In very crowded fields, such as those observed towards the 
Galactic centre, the observed objects can be the blend of 
several stars. This blending effect is a major source of con- 
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cern for the interpretation of microlensing searches. This 
is the reason, we recall, that led, in evaluating the optical 
depth, to the choice of considering only bright sources for 
which one expects the blending effects to be alleviated. 
The multiple effects of blending are supposed to roughly 
bal ance each other when ev aluating the optical depth (see 
e.g. Hamadache et al.ll2006t) . On the other hand, as blend- 
ing is expected to cause an underestimation of the eval- 
uated event duration, we may question its relevance with 
respect to our results. 

The extent to which blending contaminates the 
sample of red clunip giants is a subject of debate 
(IPopqwski et al.ll2005l: Sumi et al.|[2006l : iHamadache et al 



2006t ISmith et al.l 120071 ). For the present analysis we re- 
mark that both MACHO and EROS evaluate the opti- 
cal depth without including the effect of blending, while 
OGLE, who find this effect to be relevant within their 
data set, use blended fits. Throughout our analysis we 
have used the reported values of the duration according 
to this choice. 

IPopowski et all (|2005h . for the MACHO collaboration, 
identify an "extremely conservative" subset of events for 
which they evaluate the blend fraction to be very close to 
1. We carry out our analysis on thi s subsample, composed 



Popowski etldl (|2005l ) 



of 22 events within the CGR. As m 
we introduce an overall normalisation factor for the mi- 
crolensing rate equal to the ratio of the number of events 
in this restricted sample to that of the complete sample. 
This is coherent with the purpose of the analysis, where 
one wants to test whether blending substantially affects 
the event parameters, and the derived quantities such as 
the optical depth and the microlensing rate, while assum- 
ing that it does not change the number of detected events. 
As a result we find somewhat broader contours, because of 
the smaller number of events, with the brown dwarf slope 
unbounded at its lower tail, but otherwise fully compat- 
ible with our previous results. This is in agreement with 
our previous discussion. Indeed, the average observed du- 
ration for this sample turns out to be similar to the full 
CGR sample (^e) = 22.3 d but 5 out of the 6 very short 
duration events are excluded. 



On the other hand ISumi et al.l (|2006'). for the OGLE 



collaboration, worked the other way round. They repeated 
their analysis assuming no blending, finding a new sample 
of 48 microlensing candidates, with an average duration 
roughly 20% shorter than in the 32 events sample. Within 
this new sample there are 3 candidates with t^ < 2 d, 
and this is of course relevant in the view of our previous 
discussion. However, these candidates are strongly affected 
by blending. We prefer, therefore, not to include them in 
our analysis. Our likelihood analysis carried out on this 
45-event subsample turns is compatible with the previous 

one. 

Finally, we recall that EROS (Hamadac he et al.ll2006h 
conclude that blending docs not affect significantly their 
results, and comment the apparent discrepancy with the 
result obtained by OGLE on this issue on the basis of their 
different choice for the threshold value of the amplification 
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Fig. 8. Cumulative duration distributions (P) for (top to 
bottom) the MACHO, EROS and OGLE sample of events. 
Superimposed, the predicted distribution for the fiducial 
model. 

(in particular, OGLE consider also less amplified events 
with respect to both EROS and MACHO analyses). 

In conclusion, given the available data sets, blending, 
though relevant, does not appear to significantly affect our 
results. 



4.5. Long duration events 

In Fig. [8] we show, for the three complete sets of events 
we consider (MACHO, EROS and OGLE), the cumula- 
tive duration distribution and the expected cumulative 
distribution for the fiducial model (averaged as for the KS 
analysis in Sect. 14.275]) . evaluated at the IMF slopes that 
maximise the likelihood. As for the smaller CGR sample 
in Fig. [71 we note the rather good agreement especially 
for short duration events. For both EROS and MACHO 
data sets we also observe a systematic excess of long dura- 
tion events. This turns out to be, however, only marginally 
significant. For instance, the models predict 10% of events 
with t-E > 51 d, compared with ^ 15% of the observed 
events. 

5. Conclusions 

We have considered the sample of microlensing 
events observed towards the Galactic bulge with red 
clump gi ant so urces reported by the th e MA CHO 
[Popowski et al.l 12005). O GLE (|Sumi et all l2006l) and 
EROS ( Hamadache et al. 200(]| ) collaborations to place 
constraints on the bulge mass function. In particular, 
through a likelihood analysis, we have studied the slopes, 
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Q^BDjC^MS) of a power law mass function in the brown 
dwarf (0.01 - 0.08) Mq and main sequence (0.08 - 1.) Mq 
mass ranges. 

For our fiducial model, comparing to the CGR sample 
of 42 MACHO events, we obtain aMS = 1-7 ± 0.5. This 
result co nipares well to that o btained in the (0.15—1) Mq 



range bv iZoccali et al 



(|2000l ). a ~ 1.3. The slope in the 
brown dwarf range turns out to be less well constrained, 
Qbd = 1.6 ±1. Overall our maximum likelihood results in- 
dicate a rather significant contribution of low mass lenses, 
with ~ 30% of the events to be attributed to brown 
dwarf s. The last result is in agreement with previo us anal- 
yses (|Han fc Gouldl Il996t ICrenacher et al.lll999l ). while 
our value for the main sequence slope, somewhat smaller, 
may be explained because we are using a more suitable 
sample of red-clump-source events. 

The analyses of the EROS and OGLE data sets give us 
somewhat different results. We derive a smaller value for 
the slope in the main sequence range, although compati- 
ble with the MACHO data set result, but we obtain only 
an upper limit for the slope of the brown dwarf popula- 
tion. This behaviour finds its explanation in the different 
observed timescales. In particular, very short timescale 
events {tE < 5 d), only observed in the MACHO data 
set, are essential to constrain the brown dwarf mass func- 
tion. The lack of short timescale events has be en noticed 
in the EROS analysis (jHamadache et al. I l200d) . In all of 



the experiments, the detection efficiency of short dura- 
tions events is extremely low, rendering the analysis in 
the brown dwarf regime difhcult and, therefore, making 
the result less robust. 

More reliable constraints on the mass function may 
come from a better understanding of the bulge model, but 
especially, as already stressed , by improving the statistics 
of observed short duration events. 

Furthermore, we have carried out an analysis on the 
optical dep th. The agreement w ith the expected values is 
recognized ()Han fc Gouldll2003h . Here we have considered 
the profile of the expected optical depth as compared to 
the observed one, finding a good agreement for both the 



model s we have considered, the model E2 of [ Stanek et al 
( 1997t ) and the model G2 of lDwek et all (ll995|). Tofurthei 



constrain the bulge profile it would be useful to extend 
microlensing searches to cover a larger area in the sky 
plane. 
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Appendix A: Evaluation of the microlensing rate 

In this appendix we detail the evaluation of the differential 
microlensing rate (Sect. [3]). This quantity is directly re- 
lated to the number of expected microlensing events eval- 
uated as 



dA^,, 



NobsTobsdV 



(A.l) 



where dP is the differential rate at which a single star is 
microlensed, N^^s is the number of monitored sources and 
Tobs is the whole observation time. 

The microlensing rate expresses the number of lenses 
that pass through the microlensing tube d'^a; in the time 
interval di, for a given number density distribution ni{x) 
and velocity distribut i on fivj) of the lenses. It reads 
(|De Ruiula et al.lll99ll : lGriesdll99lh 



dr = '"^77 ^ xfivM^'vixf^^^ '-xf{vs)dHs.{A.2) 

at Is 

The two last terms account respectively for the spa- 
tial and velocity distribution of the sources, with = 
J psD^^"^ dDs- We take into account that the volume el- 
ement varies with distance as D^dDg, Dg being the dis- 
tance between observer and source, and that the fraction 
of monitored stars having a luminosity higher than a min- 
i mum detectable luni inosity, L.^,, scales as L^''^'^ ~ Dj"' 
( Kiraga fc Paczvnskil [l994i) . Throughout the analysis of 
this paper we use 7 = 0, characteristic for bright stars 
that can be considered, at least approximately, as stan- 
dard candles. 

The volume element of the microlensing tube is d^x — 
{vr^±-n)dtdS. dS = dldDi is the portion of the tube exter- 
nal surface, Di is the distance between observer and lens, 
and dl — utREda, where Re is the lens Einstein radius, 
ut is the maximum impact parameter, Vr^± is the compo- 
nent of lens velocity in the plane orthogonal to the line of 
sight (hereafter los), and n is the unit vector normal to 
the tube inner surface at the point where the microlensing 
tube is crossed by the lens. In the following 6 is the angle 
between Vr and n, with 6* G (^f i f ) (one considers only 
lenses that enter the tube). 

The velocity of the lenses entering the tube reads 



Vi = Vr 



Vt, 



(A.3) 



where Vt is the tube velocity. On the lens plane, we have 
Vt = {1 — x)Vobs + XV where a; = v^bs is the ob- 
server's (solar) velocity and Vg is the source velocity. We 
decompose both lens and source velocities into a random 
plus a bulk component v = v + Vdri ft ■ In conclusion 



Vi = Vr 



XV ^ 



(A.4) 



where we have defined the vector A so as to include all the 
bulk motion, A= {vobs - Vdrift.j) + x{vdrift.s - Vobs)- For 
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the random component we use an anisotropic Gaussian 
distribution fSect. [2?2)) . 

Looking at Eq. IA.2I we see that we start from the 
joint three-dimensional velocity and source distributions. 
However, only the distribution of the relative velocity on 
the lens plane is relevant to the microlensing rate, since 
it determines the lensing time scales via the relation 
= Re/vt- As we show below, i t is possible to ana - 



lytically evaluate this distribution ( Riffeser et al. 20061 ) 



Indeed, besides the velocity components along the los, 
both the remaining components of the source velocity can 
be analytically integrated. A final integration on the re- 
maining component of the lens velocity is not possible, 
however, as a consequence of the assumed anisotropy of 
the velocity distribution. 

The rationale of the evaluation is as follows. A first 
integration along the los for both sources and lenses leaves 
us with two Gaussian non diagonal distributions, that we 
project on the lens plane and diagonalise. This defines 
two frames on the lens plane whose axes, in general, will 
be misaligned. We then fix one of the two frames as a 
"reference" , in particular that relative to the lens velocity 
distribution, we evaluate the relative velocity distribution 
by making use of Eq. IA.4I and integrate out the source 
velocity distribution components. 

After integration along the line of sight the two- 
dimensional velocity distribution on the lens plane orthog- 
onal to the los is 



(')p.i 

^3 - 

(i)p,l Q 



(')P.2 

— 



(27r)2cr(i)p^ia-(i)p^2 



{i)p : 



(A.5) 



where the suffix [i) indicate either lenses or sources. For 
the last component, because of the projection, Vsp — x Vs 

and cTsp{i,2} = X crs,{i,2}- 

Let the principal axes of the intersection ellipse of the 
lens (source) proper velocity ellipsoid with the lens plane 
be {xi^i,xi^2} {{xs^i, Xs.2}), hereafter we refer to the for- 
mer frame as OL, lu the angle between xi^i and x^^i and 
{vspi,i, Vspi^2} the source velocity components in OL. lu, as 
well as the values for the projected dispersion velocities, 
are fixed by the geometry, once the los has been chosen (w 
varies up to ~ 8°, increasing with the Galactic latitude). 
The distribution for the relative velocity is then evaluated 
as 



P{Vr)d^Vr= J P{Vr,Vspl)(fvspld^Vr = 



f{Vspl)f{vip)5{Vr-{vip - Vspl 
1 



A))d?vipA^Vspi(i^Vr^ 



where 



_g-(«,+A)-E-(«,+A)j2^ 



y2 



Sa + Sb 

Sa - Eb 



2 I 2 I 2/2 I 2 \ 



Fig. A.l. The microlensing tube cross section and the 
angles involved. 



2\l/2 



^)Y) 

^^Ki--?2)«i-<2)sin^-)V2 



(A.7) 



Here both v^. and A must be evaluated in OL. The angle 
between A and xi^i is fixed by the geometry once we assign 
the los. Moreover, as Fig. lA. II shows. /? = 7r + a — 6'is the 
angle between Vr and A. 

Finally, after moving to polar coordinates on the lens 
plane, {vr,0), we reach the following expression for the 
microlensing rate 



dr = 2f{^,)p, 



I, 



RE{fJ-,Di,Ds)ut X 



P{vr, a) v'^dvr d/idZ?/ dDs da , 



(A. 



where we have exploited the periodicity of the trigonomet- 
ric functions involved in P{vr) to analytically integrate 
over 6 (this provides the factor '2'). Because of the as- 
sumed anisotropic velocity dispersions, the analytical in- 
tegration over a is not possible and the dependence on 
this variable survives in the relative velocity distribution. 
The expression for the differential rate dF/d^E is easily 
obtained from Eg. lA. 81 using the relation Vr = Re/Te. 

In Eq. IA.8I we have introduced the dependence of the 
lens number distribution on the mass of the lens, fj,, with 
the usual "factorisation hypothesis" stating that the lens 
mass di stribution is independe nt of the lens spatial distri- 
bution (jPe Ruiula et al. I ll99lh . pi is the lens spatial dis- 
tribution, f{fJ.) the lens mass function that we normalise 
as follows 



(A.6) 



/(^)A*dA*=^ 



(A.9) 



where po,i is the central density. 
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